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ABSTRACT 

The impact of new laboratory data for the reaction 
BrO + H0 2 ---> HOBr + 0 2 in the depletion of global 
stratospheric ozone has been estimated using a one- 
dimensional photochemical model taking into account the 
heterogeneous reaction on sulphate aerosols which converts 
N 2 O s into HN0 3 . Assuming an aerosol loading 2 times as 
large as the "background" and a reaction probability of 0.1 for 
the above heterogeneous reacuon, the 6 fold increase in the 
measured rate constant for the reaction of BrO with H0 2 
increases the computed depletion of global ozone produced by 
20 ppt of total bromine from 2.01% to 2.36%. The use of the 
higher rate constant increases the HOBr mixing ratio and 
makes the bromine partitioning and the ozone depletion very 
sensitive to the branching ratio of the potential channel 
forming HBr in the BrO + H0 2 reaction. 

1. INTRODUCTION 

Bromine compounds play a significant role in both the 
depletion of global stratospheric ozone (Wofsy et al., 1975), 
(Yung et ah, 1980) and in the perturbed chemistry which 
leads to the ozone hole formation in polar stratospheric 
regions (Mac Elroy et al.,^986). 

Assuming kinetic data available in 1980, Yung et al. 
(1980) concluded that the catalytic cycle involving the 
reaction BrO + CIO was the main cycle involving bromine 
compounds to destroy odd oxygen in the lower stratosphere. 
More specifically, the rate constant they assume for the 
reaction: 


BrO + H0 2 — >HOBr+0 2 (1) 

was too low to make the following cycle significant: 

Br + 0 3 — > BrO + 0 2 (2) 

BiO + II0 2 — > HOBr + 0 2 (1) 

HOBr + hv — > OH + Br (3) 

OH + 0 3 — > HO, + 0 2 (4) 


net: 20 3 — > 30, 


Recently, Poulet et al. (1992) have reported a rate 
constant for reaction (1) which has been measured at 
LCSR/CNRS to be 6 times higher than the preferred value 
given in the kinetic data bases (De More et al., 1990). The 
new value is k, = 3.3 10" u cm 3 . molecule' 1 . s' 1 instead of k t = 
5.0 10 12 cm 3 . molecule 1 . s' 1 . In their modelling of the influence 
of this new kinetic data on the stratospheric chemistry, they 
show that the reduction of global ozone column density 
produced by 20 ppt of bromine is increased from 1.14% to 
1.45%. They conclude that the above cycle is no longer 
insignificant to deplete ozone. Besides, they point out that a 
possible channel for reaction (1) yielding HBr and 0 3 as 
product* CQuld decrease significantly the ozone depletion 
assuming the new kinetic data. They found that with a 10% 
branching ratio for the HBr forming channel of reaction (1), 
20 ppt of bromine would deplete the ozone column by only 
0.6%. This channel may occur at low temperatures, but this 
has to be established by laboratory studies. 

The present paper describes the impact of the new 
kinetic data on the bromine chemistry of the lower 
stratosphere in the presence of aerosols. 

2. IMPACT OF THE BrO + H0 2 REACTION IN 
PRESENCE OF AEROSOLS 

A ID steady state photochemical model has been used 
to estimate the effect of the new kinetic data on both the 
bromine partitioning in the stratosphere and the global ozone 
depletion due to bromine compounds in the presence of an 
aerosol layer. This layer is assumed to reduce the 
concentration of the NO x compounds. The effect of the 
possible occurence of the HBr forming channel in reaction (1) 
is also investigated. 

The ID model has been recently described 
(Ramaroson et al., 1992). This model includes the species of 
the O x , HO x , NO y> C10 y and BrO Y families and the source 
species: N 2 0, H 2 0, CH„ CO, CFC S , CH 3 Cl and CH 3 Br. The 
vertical temperature and total concentration profiles are taken 
from the U.S. Standard Atmosphere (1976). The kinetic and 
photochemical data used are essentially those recommended 
in the last NASA-JPL report (De More et al., 1990). Those 
concerning the bromine species are reported in Table I.The 
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following heterogeneous reaction on sulfate aerosols has been 
taken into account: 

N 2 O s + H z O (sulfuric acid solution) — -> 2HN0 3 (5) 

The reaction probability has been assumed to be y - 0.1 
which is consistent with measurements of Hanson and 
Ravishankara (1991). A sulfate aerosol area profile twice as 
high as the ’’background" 1979 values over Laramie, 
Wyoming (Hofmann and Solomon, 1989) has been adopted 
in the calculations (Rodriguez et al., 1991). The stratospheric 
aerosol layer is assumed to be saturated and all the HN0 3 
produced in reaction (5) is immediately released in the gas- 
phase. 

The steady-state vertical distribution of the source 
species and families are computed by iterations from 0 to 60 
km, by step of 1 km. The diurnal averaged production and 
loss terms needed are calculated before each iteration, which 


Tabic I: Reactions involving bromine compounds and rate coefficients 


CH,Br + hv -» CH, ♦ Br 
BrO ♦ hv — ♦ Br ♦ O 
BrONO, ♦ hv -4 Br + NO, 
HOBr 4 hv -4 OH ♦ Br 
BrCI * hv -4 Br + Ci 
CH,Br ♦ OH — * products 1 * 

Br * O, — i BrO + O, 

BrO + O -* 8 r ♦ O, 

BrO ♦ CIO -+ Br + ClOO 
BrO 4 CIO -4 Br 4 OCIO 
BrO 4 00-1 BrO 4 O, 

BrO 4 NO -4 Br 4 NO, 

BrO 4 BrO -» Br 4 Br 4 O, 
Br 4 HO, -4 HBr 4 O, 

HBr 4 OH -* Br 4 H,0 
BrO 4 HO, -4 HOBr 4 O, 
BrO 4 OH -4 Br 4 110, 

HBr 4 O -4 Br 4 OH 
Br ♦ CH ,0 HBr 4 CHO 
BrO 4 NO, 4 M — » BrONO, 4 


(a) 

(a) 

(a) 

(*) 

(a) 

iJ& 10 ° exp(-450/T) cm’ molecule 1 s' 1 
1.7 10 11 exp<-aOO/T) cm* molecule 4 s' 4 

3.0 10 “ era’ molecule * s'* 

2L9 10 “ exp(H0/T) cm’ molecule' s * 
Ij6 10 “ «xp(<30rT) cm* molecule 4 s‘‘ 
5J Iff" exp{170/T) cm’ molecule^ s'* 
tS Iff" rxp(240/T) cm’ molecule * s’ 1 
1 A 10 “ cxp<l50/T) cm* molecule-* « 1 
1-5 10* 1 exp(-400/T) cm* molecule' 1 s'* 

1.1 10" cm’ molecule' 1 s * 

5/) Iff" cm* molecule'* s'* 

1 JO 1C" cm* molecule-* s'* 

SJt Iff" exp(-lS00/T) cm' molecule * j * 
1.7 Iff" exp<-at>0/T) cm* molecule' 4 »"* 
k, = 5.2 10" (T/300) 14 cm* molecule ’ **• 
k. = 9.0 10 11 (TO0C) “ cm’ molecule * s’ 


(a) : sec text for the photodis social ion rate caJculations 

(b) : considered in the model as the limiting step for the Br production from CH,Br 

(c) - k = ) 06 1 1 + Cog^MMI/k-)) 1 )•' 


requires the diurnal evolution of the short-lived species. 
Starting with realistic vertical distributions, convergence is 
obtained after 12 iterations. In the calculations, a total 
chlorine content of 3.2 ppb and a total bromine content of 20 
ppt have been used. 

Figure 1 compares the computed vertical distributions 
of N0 2 , H0 2 and CIO at noon if reaction (5) is included in 
the photochemical scheme (case b) with the same vertical 
distributions if (5) is not included (case a). Reaction (5) leads 
to an important sink of NO x (NO + NO z + N0 3 + 2 x N 2 O s ), 
as well as to an additional source of HO x due to the 
subsequent photolysis of HN0 3 . Chlorine compounds are also 
affected by this reaction. High levels of HO Xs through the 
reaction of OH with HC1, increase the concentration of the 
reactive chlorine species (Cl + CIO + C!ON0 2 ). Low values 
of N0 2 shift the partitioning between CIO and C10N0 2 in 
favour of CIO. Figure 1 shows that reaction (5) decreases the 
N0 2 concentration by 45% and increases the concentration 
of H0 2 and CIO by respectively 36%= and 1 00%, at 20 km. 



Figure l: VerLicaJ distribution, at noon, of the concentrations of NO z , H0 2 
and GO. Case (a): reaction (5) is not included. Case (b): reaction (5) is 
included. 


These results are globally consistent with the 2D 
model calculations of Rodriguez et al.( 1991 ). We may assume 
in consequence that the results of our ID calculations 
concerning the influence of the new kinetic data for reaction 
(1) on both the partitioning of bromine and the ozone 
depletion are realistic at least globally. 

Figure 2 shows the combined effect of the high rate 
constant k, and of the occurence of reaction (5), on the 
diurnal averaged concentration of bromine species. 
Concentrations of BrO and BrON0 2 are mainly affected by 
reaction (5) in the same manner than chlorine compounds. 
HOBr is affected by both the reaction (5) and the high 



Figure 2: Vertical dist:::uiiion of the diurnal averaged concentration of HOBr, 
BrONO z and BrO. Case (a): reaction (5) is not included and k 1 = 5.0 10 12 
cm 5 .molcculc l .s ’. Case (b): reaction (5) is included and k, = 3.3 10 “ 
cm 3 . molecule '.s’ 1 . 
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constant k r At 20 km its averaged concentration is increased 





by a factor 20. 

The effect of the new rate constant on the efficiency 
of the catalytic cycles involving bromine compounds to 
deplete global ozone has been also computed. Four cycles 
have been assumed: the two cycles already mentionned which 
involve BrO + CIO (cycle I) and BrO +H0 2 (cycle II)together 
with those which involve BrO + BrO (cycle III) and 
BrO + O (cycle IV). Figure 3 shows the odd oxygen 
destruction rates by cycles I to IV using the new kinetic data 
for reaction (1). The odd oxygen destruction rate by all the 
other cycles is also shown for comparison. The calculations 
show that the new rate constant has no effect on cycles I, III 



Figure 3: Odd oxygen destruction rate by catalytic cycles I, II, III and rv 
defined in the text, as a function of altitude, ALL OTHERS refers to the sum 
of the destruction rates due to the cycles which do not involved bromine 
species. Reaction (5) is included in the phochemical scheme. The destruction 
rates for cycles I, III, IV and ALL OTHERS are computed with k^ =3.3 10 " 
cnP.molecule'Vs' 1 . For cycle II, case (a): k, = 5.0 10* 12 cnP.moleculeTs \ case 
(b): k, = 3.3 10 " cm 3 . molecule 1 .s' 1 

and IV. Figure 3 shows that the rate of cycle II increases by 
a factor 6 similarly to the increase of the new rate constant. 
This was expected because reaction (1) is the limiting step of 
cycle II. 

The destruction rate by cycle II is lower than the 
destruction rate by cycle I by about a factor 4.5 at 20 km. As 
it can be seen from Figure 4 in Poulet et al.( 1 992), this factor 
was only 3 when reaction (5) was not taken into account. The 
importance of cycle II to deplete ozone is therefore reduced 
by the presence of aerosols, relatively to cycle I. It remains 
nevertheless significant. 

To quantify the impact of the new rate constant to 
deplete ozone, the percentage of ozone decrease versus 
altitude produced by 20 ppt of bromine has been calculated 
(Figure 4) using successively the low and the high value for 
the rate constant k t . From these calculations, it is found that 
the new rate constant leads to an increased reduction of the 
ozone column density from 2.01% to 2.36%. This is 
significantly higher than the reduction of 1.14% and 1.45% 


reported by Poulet et al.(1992). 

The channel of reaction (1) yielding HBr and 03 is 
not unlikely, mainly at low temperature. Laboratory studies 
have to establish the branching ratio of this channel. As 
already pointed out (Poulet et al., 1992), even a small value 
of this ratio will lead to a large increase of the concentration 
of HBr in the lower stratosphere making HBr one of the most 
abundant species of the bromine family (Figure 5). Figure 4 
shows also, the percentage of ozone decrease versus altitude 
produced by 20 ppt of bromine for a 10% branching ratio of 
the HBr forming channel and the higher value of We 



Figure 4: Percentage of ozone depletion produced by 20 ppt of bromine, as 
a function of altitude. Heterogeneous reaction (5) is included. Case (a): k x = 
5.0 10‘ 12 cm\s' 1 , branching ratio for the Hbr forming ratio is zero. Case (b): 
k, = 3.3 10 " cm 3 .s' 1 , branching ratio is zero. Case (c): k, = 3.3 10 " cm 3 .s' 1 , 
branching ratio is 10%. 

observe a strong decrease of the ozone reduction. The 
calculations lead to a decrease of the ozone column density 



Figure 5: Vertical distribution of the concentration of bromine species, at 
noon. Reaction (5) is included, k,= 3.3 10 " cm 3 . molecule' 1 s' 1 , the branching 
ratio for die HBr forming channel is 10%. 
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of only 0.75%. These results were expected because HBr is 
a more efficient reservoir of active bromine than HOBr. We 
can observe that the inclusion of reaction (5) in the 
photochemical scheme increases the difference in the 
reduction of the computed ozone column assuming 10% and 
0% for the branching ratio of the HBr forming channel. Using 
reaction (5) and the higher value of k lt the decrease of the 
reduction of the ozone column is from 2.36% to 0.75%. 
Neglecting reaction (5), the decrease was only from 1.45% to 
0.6%. 

3. CONCLUSION 

The modelling results presented in this paper, 
compared with those of Poulet et al. (1992) which did not 
take into account the heterogeneous reaction: 

N 2 0 5 + H 2 0 -> 2 HN0 3 

in the aerosol layer, show that the importance of the catalytic 
cycle involving the reaction BrO + H0 2 is reduced compared 
to the cycle involving the reaction BrO + CIO, but remains 
significant. In the present calculations, the ozone column 
density reduction produced by 20 ppt of total bromine 
increases from 2.01% (with the lower rate constant kj) to 
2.36% (with the higher rate constant k t ) assuming an aerosol 
loading twice as high as the "background 1 '. 

An increase in the aerosol loading would not changed 
significantly these conclusions. Others calculations have been 
repeated assuming the surface area profile available at 
Laramie at the maximum of the El Chichon eruption 
(Hofmann and Solomon, 1989) which was about 25 times as 
large as the "background”. As already pointed out by 
Rodriguez et al. (1991) the reduction of the ozone column 
density is not too much affected by this large aerosol loading 
because the concentration of the NO x compounds relatively 
to the concentration of HNO, becomes independent of the 
above heterogeneous reaction when aerosol loading increases. 
In the presence of this large surface area profile, we calculate 
a reduction of the ozone column density of 2.6% produced by 
20 ppt of bromine with the lower rate, constant for the 
reaction of BrO with H0 2 and a reduction of 3% with the 
higher rate constant. The cycle involving the reaction BrO + 
H0 2 would be therefore significant to deplete ozone even in 
presence of a large aerosol loading if the new rate constant is 
taken into account. 

The calculations made assuming a 10% branching ratio 
for the HBr forming channel of reaction BrO + H0 2 shows an 
increasing importance of this channel if the above 
heterogeneous reaction in the aerosol layer is included and 
when the aerosol loading is increased. This channel would 
decrease the reduction of the density ozone column produced 
by bromine species. When the heterogeneous reaction was not 
included in the photochemical scheme a 10% branching ratio 
decreases the reduction of the ozone column from 1.45% to 
0.6% using the new rate constant, assuming 20 ppt of total 
bromine. When the heterogeneous reaction is used, the 


decrease is from 2.36% to 0.75% with a aerosol loading 2 
times as large as the "background” and 3% to 0.85% with a 
loading 25 times as large as the "background". 

Ozone depletion by bromine species is therefore very 
dependent on the branching ratio for the HBr forming channel 
of reaction (1). Laboratory investigations are therefore 
needed, at low temperature, to measure this branching ratio as 
well as in-situ measurements of HBr. Park et al.(1989) report 
a far-infrared measurement of 20 ppt at 28 km while Traub 
et al.( 1992) report three measurements and conclude to an 
upper limit of 4 ppt at 32 km. After Traub et al.(1992), Park 
and his co-workers now agree that 20 ppt is in fact a tentative 
upper limit. Our calculations (Figure 5) give 6 ppt at 32 km, 
assuming a 10% branching ratio. A few percent for this 
branching ratio would not be therefore inconsistent with the 
measurements of Traub et al.(1992). Other measurements are 
obviously needed. 
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